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Introduction 

The  binary  gas  mixture  CF4/02  is  widely  used  in  industrial  production  of  semiconductor 
devices  to  increase  the  etching  rate  of  silicon  wafers.  The  number  of  papers  devoted  to  mathe¬ 
matical  modeling  of  silicon  etching  in  CF4/02  plasma  is  relatively  small.  Moreover,  the  numeri¬ 
cal  modeling  of  silicon  etching  with  simplifying  kinetics  of  the  gas  phase  and  heterogeneous 
chemical  reactions  predicts  a  maximum  of  the  etching  rate  near  the  50%  fraction  of  02  in  the 
parent  gas  mixture  [1].  This  is  inconsistent  with  experimental  data  where  the  maximum  etching 
rate  detected  was  at  the  15-30%  range  of  02. 

In  this  paper,  the  results  of  numerical  optimization  of  a  radial  flow  plasma-chemical  reac¬ 
tor  depending  on  the  binary  gas  composition  CF4/02  are  presented.  Quantitative  estimations  of 
oxygen  chemisorption  and  adsorption  of  CF2,  CF3  contributuons  to  the  etching  rate  are  obtained. 

Mathematical  model  of  radial  flow  plasma-chemical  reactor 

The  authenticity  of  received  numerical  results  depends  on  the  choice  of  chemical  kinetics 
of  plasma  etching.  Using  comparison  with  experimental  data  for  plasma  CF4/O2  in  [2]  the  mini¬ 
mum  set  of  chemical  reactions  was  derived  which  satisfactorily  predicted  the  concentrations  of 
stable  reaction  products  COF2,  C02,  CO  and  F.  In  [3]  this  basic  set  of  reactions  was  added  by 
heterogeneous  reactions  of  species  F,  CF2,  CF3,  CF4,  O  to  describe  the  adsorption  of  CF2,  CF3. 
The  improved  chemical  kinetic  model  [3]  including  12  components  -  F,  F2,  CF2,  CF3,  CF4,  C2F6, 
O,  02,  CO,  C02,  COF,  COF2  was  used  in  present  calculations.  The  heat-mass  transfer  in  the 
reactor  was  calculated  on  the  base  of  two-dimentional  mathematical  model  of  nonisothermal 
reactor  [4]  with  taking  into  account  the  chosen  multicomponent  chemical  kinetics.  The  heat  ra¬ 
diation  in  the  gas  mixture  was  determined  in  optical  thin  layer  approximation.  To  evaluate  the 
emissivity  of  gas  mixture  the  exponential  model  of  band  adsorption  [5]  describing  the  spectral 
adsorption  of  multiatomic  molecule  was  employed.  The  experimental  data  shown  that  the  radi¬ 
cals  CF2,  CF3,  C2F4,  C2F5  adsorbed  on  silicon  surface  form  the  recombination  cites  acting  as  a 
barrier  preventing  an  etching.  The  quantitative  estimations  for  silicon  etching  in  pure  CF4  pre¬ 
sented  in  [4]  have  shown  that  the  effect  of  CF3  adsorption  on  etching  rate  does  not  exceed  5%. 
In  CFVCL  plasma  the  oxygen  reacts  with  CFX  radicals  preventing  their  recombination  with  F  and 
reducing  their  concentrations.  Nevertheless  it  is  of  interest  to  estimate  the  influence  of  adsorp¬ 
tion  layers  on  the  etching  rate  also  in  this  case.  It  is  of  great  significance  during  silicon  etching  a 
formation  of  active  centers  of  fluorine  and  oxygen  atoms  on  silicon  surface.  This  process  was 
described  in  approach  represented  in  [6].  The  governing  equations  were  numerically  solved  on 
the  base  of  iterative  finit-difference  splitting-up  schemes  with  stabilizing  correction. 

Results  and  discussion 

The  calculations  of  loaded  reactor  have  been  completed  for  several  values  of  gas  flow  rate 
Q  and  wafer  temperature  Ts  presented  in  Table  1.  The  pressure  in  etching  chamber  of  reactor 
was  equal  to  p  =  0.5  torr.  The  average  electron  density  was  essumed  equal  to  ne  =  6  x  109,  cm3. 
The  temperature  of  reactor  walls  -  Tw  =  300  K.  The  02  persentage  fraction  in  CF4/02  feed  gas 
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Fig.  1.  Isolines  of  stream  function  in  the  radial  flow  reactor.  The  arrowed  line  shows  the  directions  of 
gas  flow  in  the  reactor:  a  -  “inflow”,  b  -  “outflow”.  Processing  regime:  p  =  0.5  torr,  Q  =  800  cm3/min, 

Ts  =  300  K,  30%  fraction  of  02  in  CF4/O2  input  gas  mixture. 

mixture  varied  in  the  range  10-90%.  In  order  to  compare  the  results  of  our  calculations  with  the 
results  [3]  the  constructive  dimensions  and  operating  regimes  were  taken  as  in  [3].  The  wafer 
temperatures  coincided  with  data  [7].  Two  gas  flow  directions  were  examined:  “inflow”,  when 
the  gas  flow  was  directed  to  the  center  of  reactor,  and  “outflow”  -  with  reversed  direction. 

Flow  structure  and  heat  transfer.  The  isolines  of  stream  function  in  radial  flow  reactor 
are  presented  on  Figure  1.  The  direction  of  feed  gas  flow  determined  a  peculiarities  of  flow 
structure  in  the  reactor.  In  the  “inflow”  case  the  flow  structure  in  the  etching  chamber  has  no 
large  vortices.  The  velocity  distribution  in  the  reactor  is  weakly  dependent  on  gas  flow  rate  and 
wafer  temperature.  Under  “outflow”  direction  a  vigorous  gas  stream  is  injected  to  the  etching 
chamber  and  then  is  spreaded  out  along  the  RF-electrodes.  Near  the  inlet  tube  arranged  in  the 
centre  of  lower  electrode  the  recirculating  vortical  zone  is  formed.  The  radial  dimension  of  vor¬ 
tical  zone  is  increased  with  rising  of  gas  flow  rate  or  decreasing  of  02  percentage  fraction  in 
mixture.  The  ratio  of  characteristic  velocities  w,  to  w0  at  inner  and  outer  reactor  edges  approxi¬ 
mately  consists  of  w,- /  Wo  ~  102. 

For  both  feed  gas  flows  the  temperature  of  gas  mixture  monotonically  rises  along  the  flow 
direction.  The  maximum  gas  temperature  localizes  near  the  wafer.  For  “inflow”  direction  it  lo- 
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Fig.  2.  Isolines  of  fluorine  concentration,  Cp  x  10"8  Mol/cm3.  Processing  regime:  p  =  0.5  torr, 
Q=  200  cm3/min,  T,  =  300  K,  40%  fraction  of  O?  in  CF4/O2  input  mixture. 
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Fig.  3.  Distribution  of  full  flow  of  active  species.  Processing  regime:  see  Fig.2. 

cates  in  the  range  r  =  2.9-11.8  cm  (Fig.  l,a).  In  opposit  case  it  correspondingly  falls  within  the 
range  r  =  15.0-21.4  cm  (Fig.  l,b).  The  typical  differences  of  gas  temperature  in  the  reactor  are 
found  in  interval  42-177  K.  The  rise  of  CF  fraction  of  in  the  range  10-90%  in  gas  mixture  re¬ 
duces  the  gas  temperature  on  9-40  K,  that  is  explained  by  the  decreasing  of  emissivity  of  gas 
mixture.  With  the  rising  of  wafer  temperature  the  radial  temperature  gradients  are  increased  that 
brings  down  an  etching  wafer  unuformity.  The  distribution  of  full  heat  flow  for  the  most  part  of 
reactor  coincides  with  the  direction  of  convective  transfer.  Over  the  wafer  the  axial  component 
in  the  full  heat  flow  is  defined  by  the  heat  conduction  and  heat  radiation  flows.  The  contribu¬ 
tions  of  heat  conduction  and  heat  radiation  flows  in  the  hill  heat  flow  are  relatively  small. 

Mass  transfer.  The  distribution  of  fluorine  concentration  in  the  reactor  depends  on  many 
factors:  operating  regime  of  reactor,  gas  composition,  direction  of  feed  gas  flow  and  etc.  Figure 
2  illustrates  typical  distribution  of  fluorine  concentration  in  the  reactor.  If  the  feed  gas  flows  to 
the  center  of  reactor  the  fluorine  concentration  monotonically  rises  along  the  flow  direction  ow¬ 
ing  to  the  convective  transfer  (Fig.  2, a).  The  maximum  of  concentration  is  located  near  the  up¬ 
per  electrode  at  r  =  0-2.1  cm.  If  the  feed  gas  flows  from  a  center  of  reactor  the  distribution  of 
fluorine  concentration  is  also  monotone  (Fig.  2,b).  Flowever  at  any  fixed  cross  section  by  the 
height  of  reactor  the  fluorine  concentration  has  the  local  minimum  in  the  range  r  =  16.9-20.0 
cm.  The  position  of  concentration  maximum  is  removed  along  the  upper  electrode  with  variation 
in  the  operating  regime.  Figure  3  illustrates  the  distribution  of  full  flow  of  active  species.  The 
direction  of  fluorine  mass  transfer  in  middle  part  of  reactor  height,  where  |v|  »  w0,  coincides 
with  the  direction  of  gas  flow  (v  -  velocity  vector).  Near  the  RF-electrodes,  where  the  velocity 
of  gas  flow  is  small,  the  fluorine  mass  transfer  is  realized  by  a  concentration  diffusion.  Owing  to 
intensive  forced  convection  the  nonuniform  distribution  of  full  flow  of  active  species  on  the 
wafer  is  arised.  Because  of  that  the  full  flow  of  active  species  deviates  from  the  normal  to  the 
wafer.  Finally  as  one  can  see  from  Fig.  3,  the  mass  transfer  of  fluorine  atoms  to  the  wafer  sur¬ 
face  depends  strongly  on  the  convective  transfer  in  the  reactor. 

Etching  rate  and  uniformity.  The  local  spontaneous  etching  rate,  its  average  value  and 
uniformity  idex  were  calculated  by  the  formulas  [7,  8] : 


Vs 


Vs  =  1.81- 1010(1-  90  -  9Cf2  ~  &cf3  )ksxrC, , 
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where  Vs,  Vs/,  Vs2  -  local,  minimum  and  maximum  spontaneous  etching  rates;  9q  -  fraction  of 
silicon  surface  covered  by  oxygen  atoms,  9q\:2  >  $cf3  “  fractions  of  silicon  surface  covered  by 
adsorbed  radicals  CF2  and  CF3  correspondingly,  ks  -  etching  rate  constant.  To  simplify  a  rep¬ 
resentation  of  i9()  it  is  assumed  that  CF  »  Cq p2  ,  Cq-, .  Thus  9q  can  be  represened  as 


90 


vpap 

vGa0  +  vFaF  ’ 


where  vF,v'0  -  molar  flux  of  fluorine  and  oxygen  atoms  to  the  wafer;  aF  =  0.7,  a0  =  0.2  - 
sticking  coefficients  of  oxygen  and  fluorine  atoms  on  silicon  [6],  and  the  denominator  shows  a 
competition  of  chemisorption  of  oxygen  and  fluorine  atoms  on  silicon.  The  explicit  expressions 
for  9q\;2  ,  i9(-f  are  not  presented  here  because  of  its  awkwardness.  They  are  similar  to  used  in 

[3],  but  here  in  addition  we  take  into  account  an  oxygen  chemisorption. 

The  average  etching  rate  dependence  on  percentage  fraction  of  02  in  parent  gas  mixture 
has  the  maximum  near  40%  of  02  in  the  range  of  gas  flow  rate  100-400  cmVmin  (Fig.  4,  curve 
1).  With  increasing  of  gas  flow  rate  up  to  800  cmVmin  the  maximum  of  etching  rate  moves  to 
30%  of  02  (Fig.  4,  curve  2).  The  location  of  extremum  is  independent  of  the  feed  gas  flow  direc¬ 
tion  and  the  gas  temperature. 

With  increasing  02  percentage  in  CF4/02  gas  mixture  the  fluorine  concentration  rises  and 
has  the  maximum  at  40%  02.  Decreasing  CF2,  CF3,  CF4,  COF  concentrations  accompanies  by 
the  growth  of  fluorine  concentration.  The  substances  F2,  CO,  C02,  COF2  are  formed  as  the 
products  of  reactions.  It  corresponds  to  one  of  channels  of  etching  intensification  by  oxygen 
addition,  which  consists  in  replacement  of  fluorine  by  oxygen  in  fluorine -containing  radicals. 
With  oxygen  addition  more  than  40%  the  fluorine  concentration  falls  owing  to  decreasing  total 
number  of  CF2,  CF3  radicals  (depletion  of  mixture).  Thus  the  calculation  shows  that  the  pres¬ 
ence  of  etching  rate  maximum  connects  with  chemical  reactions  of  oxygen  atoms  with  CF2,  CF3 
radicals  where  the  additional  fluorine  atoms  set  free.  The  maxima  of  etching  rate  and  fluorine 
concentration  coincide  with  respect  to  02  addition  in  CF4/02  gas  mixture. 

The  etching  rate  dependence  on  fluorine  concentration  with  different  fractions  of  02  in 
CF4/02  mixture  obtained  in  calculations  has  a  hysteresis  character  (Fig.  5,  curve  1).  Because  of 
the  passivation  of  wafer  surface  by  oxygen  atoms  (increasing  i9q  )  in  the  range  of  O?  concentra¬ 
tion  more  then  40%  the  silicon  etching  rate  falls  more  quickly  then  fluorine  concentration  (see 
formula  (1)).  At  the  same  time  the  hysteresis  is  not  observed  under  the  absence  of  oxygen 
chemisorption.  It  is  illustrated  by  the  linear  dependence  of  Vs  on  CF  (curve  2  at  Fig.  5)  ob¬ 
tained  for  9q  =  0  and  any  fraction  of  02  in 
mixture.  This  fact  corroborates  with  the 
results  [9]  where  the  hysteresis  was  ex¬ 
perimentally  detected  and  the  assumption 
about  its  connection  with  oxygen  chemi¬ 
sorption  have  been  made.  The  influence  of 


Fig.  4.  Average  etching  rate  as  a  function  of  per¬ 
centage  fraction  of  02  in  CF4/02  input  mixture. 

1  —Q  —  200  cm3/min,  2  —  Q  =  800  cmVmin. 
Processing  regime:  p  =  0.5  torr,  Ts  =  300  K. 
“Inflow”  direction  was  examined.  Numerical 
results  were  marked  by  points,  solid  lines  - 
spline  approximation. 
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CF2,  CF3  adsorption  layers  on  the  etching 
rate  beginning  from  5%  02  concentration 
composes  less  1%  of  nominal  value. 

The  high  etching  uniformity  may  be 
achieved  by  controlling  the  processing  re¬ 
gimes  of  reactor  (gas  flow  rate,  tempera¬ 
ture,  gas  composition  and  etc.).  As  in  [3]  it 
is  found  that  the  minimum  etching  nonuni¬ 
formity  is  realized  at  the  small  gas  flow 
rates  (Table  1,  gas  flow  rate  100  cm3/min). 
Flere  it  is  possible  the  effective  control  of 
etching  rate  by  changing  the  ratio  CF4/02. 
The  rise  of  gas  flow  rate  causes  the  increase 
of  etching  nonuniformity  independently  of 
feed  gas  flow  direction.  Note  that  at  the 
gas  flow  rates  larger  then  200  cmVmin  the 
quality  of  silicon  etching  does  not  corre¬ 
spond  to  requirements  for  semiconductor 
wafers.  Table  1  shows  that  the  “inflow” 


Fig.  5.  Average  etching  rate  as  a  function  of  average 
fluorine  concentration  on  silicon  with  different  fractions 
of  02  in  CF4/02  mixture.  1  -  ^  0;  2  -  =  0. 

Processing  regimes:  p  =  0.5  torr,  Q  =  200  cm3/min, 
Ts  =  300  K.  “Inflow”  direction  was  examined. 


Table  1.  Uniformity  index  as  a  function  of  gas  flow  rate  and  wafer  temperature.  30%  fraction  of  O2  in  CF4/O2 
input  mixture.  “Inflow”  and  “outflow”  directions  were  examined. 


Q,  cnr’/min 

T„  K 

(Inflow) 

(Outflow) 

100 

300 

0.05928 

0.05385 

200 

300 

0.09776 

0.29991 

200 

373 

0.07574 

0.26711 

200 

473 

0.05121 

0.26986 

200 

573 

0.03737 

0.27081 

400 

300 

0.20619 

0.46885 

800 

300 

0.39062 

0.58001 

direction  is  more  preferable  because  of  uniformity  index  decreases  with  rising  of  gas  tempera¬ 
ture.  The  dependence  of  etching  nonuniformity  on  percentage  fraction  of  02  has  not  pronounced 
tendency.  But  its  high  value  at  percentage  fraction  of  02  more  then  80%  is  common  for  all  oper¬ 
ating  regimes.  At  this  feed  gas  composition  even  if  the  gas  flow  rate  is  100  cm3/min  the  large 
gradients  of  fluorine  concentration  are  appeared.  The  influence  of  oxygen  chemisorption  on 
silicon  surface  on  uniformity  index  is  weak. 

Conclusions 

The  obtained  results  allow  to  make  the  following  conclusions. 

1 .  The  oxygen  addition  allows  to  increase  the  etching  rate  in  several  times  and  it  is  an  ef¬ 
fective  factor  for  controlling  the  processig  regime.  The  appearing  etching  rate  maximum  con¬ 
nects  with  chemical  reactions  of  oxygen  atoms  with  CF2,  CF3  radicals  in  which  the  additional 
fluorine  atoms  set  free. 
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2.  The  fluorine  and  oxygen  chemisorption  on  silicon  leads  to  hysteresis  on  the  diagram  of 
etching  rate  with  respect  to  fluorine  concentration.  Under  the  absence  of  oxygen  chemisorption 
the  hysteresis  is  absent  too.  The  influence  of  CF2,  CF3  adsorption  layers  on  the  etching  rate 
composes  less  1%  of  nominal  value. 

3.  The  fluorine  transfer  in  the  radial  flow  reactor  to  the  wafer  essentially  depends  on  con¬ 
vective  transfer.  The  minimum  etching  nonuniformity  is  realized  at  the  small  gas  flow  rates  (100 
cm3/min)  where  fluorine  transfer  is  determined  by  a  concentration  diffusion  to  the  wafer.  Flere  it 
is  possible  the  effective  control  of  etching  rate  by  changing  the  ratio  CF4/O2.  The  “inflow  “  di¬ 
rection  of  feed  gas  mixture  provides  the  best  etching  uniformity  and  it  is  more  preferable.  At  the 
gas  flow  rates  larger  then  200  cm3/min  the  quality  of  silicon  etching  does  not  correspond  to  re¬ 
quirements  for  semiconductor  wafers. 

This  work  was  supported  by  the  Russian  Fundation  of  Basic  Research  (grant 
No.  01-01-00827). 
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